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I. INTRODUCTION 



An archetypical example of a QCD (hadronic) observable is the pion form factor, which is 
typical for a hard-scattering process obeying a factorization theorem [H, Hj. Consequently, at 
asymptotically large it can be cast in terms of a scale-dependent pion distribution ampli- 
tude (DA) [sl of leading twist two '^t,{x, Q^) convoluted with the hard-scattering amplitude 
of the process which contains the large external scale Q^: 



^pert^g2^) _ 87ra,(Q2)/2 ^ ^ ^^^^^ 



r_^{Q')V with /:i(q2^ 



dx . 



The nonperturbative input — the pion DA /x^) — can be expressed as an expansion over 
Gegenbauer polynomials 



1 + 5]]a2n(/U^ 



n>l 



where the asymptotic pion DA has the form 



^^(x) = 6x{l-x) 



(2) 



(3) 



while the scale dependence of coefficients a2n(/^^) is controlled by the ERBL evolution equa- 
tion 

At the one- loop level and at asymptotically large Q^, the pion form factor simplifies to 
= 8'7ias{Q^)f^/Q'^ The onset of the asymptotic regime cannot be determined 
precisely; estimates l5| show that this transition scale is of the order of 100 GeV^. 

On the other hand, at intermediate momentum transfers 20 GeV^ > > 1 GeV^, 
the situation is more complicated because of the interplay of perturbative and nonpertur- 
bative effects. The latter effects are contained in a non-factorizable part — called the soft 
contribution — so that one has to take it into account using some nonperturbative concepts, 
e.g., the method of QCD sum rules (SR) jl, 0, 11], the local quark-hadron duality (LD) ap- 
proach [g], 0], and others. Note in this context that, describing the pion form factor within 
the three-point QCD SR approach [gI, 0], the shape of the pion DA becomes irrelevant. This 
considerably reduces the inherent theoretical uncertainty of the method. The same applies 
to the LD approach, but the latter contains an additional uncertainty related to the So(Q^) 
setting for intermediate and large values of — see for a discussion in [sj. 

However, the standard QCD SR [1,01 are plagued by instabilities arising at ^ 3 GeV^, 
which are induced by those terms in the operator product expansion that are either con- 
stant or grow linearly with (see Tab. [T]). Such terms do not represent a nonperturbative 
contribution correctly. The corresponding diagrams result from the substitution of some 
propagators by constant factors that denote condensates lacking a Q^-dependence, viz., 
{T{q{z)q{0))) (g(O)g(O)). The scale dependence is retrieved by including in the cal- 
culation the contributions stemming from higher- dimension operators, like (g(O)D^g(O)), 
(g(0)(D^)^g(0)) etc., that are entailed by the Taylor expansion of the original nonlocal con- 
densate (NLC), {q{0)q{z)), being the nonperturbative part of the quark propagator. In order 
to obtain the correct large-Q^ behavior and ensure that the total condensate contribution 
decreases for large Q^, one has to resum all terms of the standard OPE bearing terms of the 
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sort (Q^/M^)". This is a rather tedious task and, therefore, we refrain from using the orig- 
inal Taylor expansion in our analysis and take instead recourse to a modified diagrammatic 
technique which makes use of new lines and vertices associated with NLC (details can be 
found in [i]). 

TABLE L Q^-behavior of the nonperturbative contribution in different QCD SR approaches. Here, 
ci, C2, C3, C4 are dimensionless constants (not depending on Q^). The abbreviations used are: LD 
for local duality, LO for leading order, and NLO for next-to-leading order, while and M"^ denote 
the vacuum quark virtuality and the Borel parameter, respectively. 



Approach 


Accuracy 


Condensates 


Q^-behavior of $ope 


Standard QCD SR [6, 7] 


LO 


Local 


ci + Q^M^ 


QCD SR with NLCs [8] 


LO 


Local -|- Nonlocal 


{c2 + QVM^) ^e-'^Q'^V^' + c) 


LD SR(M2 ^ 00) [5] 


NLO 







This paper 


NLO 


Nonlocal 





An earlier attempt to generalize the QCD SR [8|] approach by employing such NLC 
contributions turned out to be incomplete, because it was found to contain contributions 
originating from local condensates. This is related to the fact that a specific model (fT5|) for 
the 3-point quark-gluon-antiquark NLC was used in which the NLC Mi{x'^, y'^, {x — yY) are 
nonlocal only with respect to one single separation, say, x^, out of the three possible inter- 
parton separations x^, ?/^, and {x — y)"^. As a result, also this type of approach suffers from 
the same shortcomings as the standard QCD SR. In contrast, LD SR have no condensate 
contribution due to the 00 limit. The only trace of all contributing condensate 

contributions is embodied in the parameter sq, which can be derived from the LD sum rule 
for /tt. Due to the Ward identity, these sum rules are connected only at = 0, so that 
the applicability of this method to the pion form factor is actually confined to low momenta 
around <^ sq. The definition of sq at large is not settled in this approach [5|. 

In this presentation, we report upon an investigation of the electromagnetic pion form 
factor which employs QCD SR with NLC [s], |9|. This enables us to enlarge the region of 
applicability of the QCD SR to momenta as high as 10 GeV^. Moreover, we use a spectral 
density which includes terms of 0{as)- The influence of this NLO contribution to the pion 
form factor reaches the level of 20%. The remainder of this report is organized as follows. 
The next section contains the necessary ingredients of the QCD SR approach with NLC. The 
second part of this section contains also our results. Our conclusions are given in section UTTl 
where we further discuss our flndings in comparison with the available experimental data, 
lattice simulations, and other theoretical calculations. 



II. PION FORM FACTOR FROM QCD SUM RULES WITH NONLOCAL CON- 
DENSATES 



The nonlocality of the QCD vacuum, suggested in 10, 11, H, 13], is crucial for a correct 
determination of the pion DA and the computation of the pion form factor j^, 0] . For that 
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reason, let us first recall the main elements of this approach and discuss its application to 
the calculation of 3-point correlators in QCD. 



For the scalar and vector condensates, we employ the same Gaussian model as in [12, [13 



I.e., 

(g(O)g(^)) = m e-'^'l^^/« ; (g(0)7,g(^)) = ^ A, e'l^^'^^/^ , (4) 

where Aq = 2a sTt {qqy / 81. Note that above, a Fock-Schwinger string is attached in-between 
the quark-antiquark fields in order to preserve gauge invariance. But adopting the fixed-point 
(Fock-Schwinger) gauge z^A^{z) = each string reduces to unity, provided the integration 
path in the exponent is a straight line going from to z. The nonlocality parameter = (fc^) 
provides a useful measure of the average momentum of quarks in the QCD vacuum. It 
has been estimated in QCD SR l3| and on the lattice with a value around 



= 0.45 ±0.1 GeV . To parameterize the vector (V) and the axial- vector (A) quark-gluon- 
antiquark condensate, we use the expressions derived in lo| : 



{qiOh^^{-9Auiy))qix)) = {y^^x^ - g^^{y ■ x))Mi{x^,y^, {y - xf) 

+ (yf^Vu - g^,yy'^)M2{x^, ^/^ {y - xf) , 
^q{^)lhlp.{-gAy{y))q{x)) = ie^,yy^M^{x^, {y - xf) 

with 

oo 

M,{x\y\z^) = A, [[[dad(3djMa,(3,^)e(^'"-'^''-'^^")/\ (5) 







where the following abbreviation Ai^2,3 = x ("f'^, |) was used. The minimal Gaussian 
model of the nonlocal QCD vacuum is introduced by the following ansatz 

fr («, /5, 7) = 5 (a - A) 5 (/3 - A) 5 (7 - A) (6) 

with A = Ag/2. This model violates the QCD equations of motion, while at the same time 
the 2-point correlator of the vector currents is not transverse. To restore the QCD equations 
of motion and to minimize the non-transversity of the VV correlator, an improved model of 
the QCD vacuum was proposed [181] : 

(«, P, 7) = (1 + X,d, + Wy + W,) S {a - xA) S iP - yA) S (7 - ^A) . (7) 

Here A = |Ag and z = y, whereas 

Xi = +0.082 ; X2 = -1.298 ; X3 = +1.775 ; x = 0.788 , (8a) 
Yi = -2.243; ^2 = -0.239; 1^3 = -3.166; ?/ = 0.212 . (8b) 

These parameters satisfy the supplementary conditions 

12(X2 + F2)-9(Xi + Fi) = l, x + y = l, (9) 

following from the QCD equations of motion. 
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The Borel SR for the pion form factor, based on the 3-point AAV correlator, was consid- 
ered for local condensates in jl, 0|, whereas the NLC case was treated in [8]. This way, one 
obtains the following SR 



so so 



dsi ds2 Psisi, S2, Q'^) e 



-(si+S2)/M2 



+$G(g^M2)+$^,-,)(Q^M2). (10) 







Note that as long as the condensate terms and ^{qq) are not specified, this SR may have 
a local or nonlocal content. On the other hand, the perturbative 3-point spectral density is 
given by 



(5l,^2,Q') + 



An 



Ap(^)(si,S2,Q2) 



(11) 



The leading-order spectral density P3°''(si, S2, t) is known since the early eighties jl, 0]- As 
regards the next-to-leading order (NLO) spectral density Ap3^''(si, S2, Q^), it has been ob- 
tained quite recently [l9|. The phenomenological side of the SR contains the contribution 
which stems from higher resonances, modeled via the spectral density 

Phr(Si, S2) = [1 - d{si < So)9{s2 < So)] PsiSi, S2, Q^) (12) 

and using the continuum-threshold parameter sq. In order to irnprove the low-scale behavior 
of the pion form factor, we adopt a scheme, developed in 0, 23], and employ an analytic 



running coupling [21 



asm = ^ 



^^QCD 



bo Vln(gVA^cD 

with bo = 9 and Aqcd = 300 MeV. 

r discussion to follow, we 
local-condensate case the following expressions 



^^QCD 



For our discussion to follow, we use for the nonperturbative terms $g and ^{qq) 

0,9: 



$|5'=(m2) 



{asGG) 
12 7rM2 



$;->(g2,M2) 



104 A 







M4 



1 + 



2Q' 
13 M2 



(13) 



in the 



(14) 



These expressions, that are used in the standard QCD SR for the pion form factor, show 
a wrong behavior at large Q^: (i) The quark contribution contains both a linearly growing 
term as well as a constant one. (ii) The gluon contribution is just a constant. On the 
other hand, the perturbative term on the right-hand side of Eq. (fTOl) behaves at large 
like so/Q^ or M^/Q^. Hence, the SR becomes unstable for > 3 GeV^. But using the 
generalized QCD SR with NLC, [8|, this deficiency should not appear. Alas, even this 
approach has a dark side, because it still uses in the analysis of the pion form factor a 
partially local parameterization of the quark-gluon-antiquark NLC. Indeed, the following 
parametric functions (jSD have been used in [s] (A = A^/2): 



/f^(«,/3,7) 



X. 



S{a- XiiA) 5{P - Xi2A) 5 (7 - Xi3A) 

0.4 0.4 
1 0.4 
0.4 0.4 



(15) 
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The absence of nonlocality effects in ([5]) for the quark- ant iquark separation y"^ {i = 1) and 
also for the (anti)quark-gluon separations and (x — y)"^ {i = 2, 3) is revealed by the zero 
elements in the matrix Xij. 

Note that the NLC contributions to the pion form factor, entering the SR (ITU]) , can 
still be used in connection with an improved version of the quark-gluon NLC because the 
expressions obtained in [sl have the form of a convolution in the a-representation of the 
NLC distribution functions 7) with model-independent coefficient functions. In the 

present work we apply the minimal ([6]) and the improved ([7]) Gaussian models of NLC. The 
contribution from the vector condensate to reads 



AMM^g^) = ^ f2 + -J^)exp 



M4 V 2M2-A2 



2 M2 (2 M2 - A2) 



(16) 



This term indeed vanishes for large and is controlled by the nonlocality parameter A^. 
The larger A^, the faster this contribution decreases with Q^. The explicit expressions for 
the other condensate contributions are omitted here, but their schematic dependence can 
be read off from Table [B 

The pion form factor FTr{M'^,so), as a function of two additional parameters (Borel 
parameter) and sq (continuum threshold), is given at each fixed value of by SR Eq. (fTOl) . 
The parameter sq marks the boundary between the pion state and higher resonances {Ai, 
it', etc.). We select its value at each momentum transfer by applying to the function 
Ft^{M'^,so) the minimal-sensitivity condition with respect to the auxiliary parameter in 
the fiducial interval of the SR. These intervals and the value of the pion decay constant 
for each considered NLC model, notably the minimal and the improved Gaussian one, are 
taken from the corresponding 2-point NLC QCD SR^^ = 0.137 GeV^, Ml = 1 GeV^ and 



= 1.7 GeV^ for the minimal Gaussian model 13|, whereas for the improved one 18 
one has = 0.142 GeV^ Ml = I GeV^ and Ml = 1.9 GeV^ The continuum threshold 
Sq^((5^), which minimizes the dependence of the right-hand side of flTOj) . is fixed by the 
root-mean-square deviation x'^iQ'^^ ^o), Eq. flA.ip . in the Borel-parameter interval M^ e 
[Ml,Ml] at each value of Q^. The results of this procedure are shown in the left panel of 
Fig. [H Both models generate approximately constant values of so(Q^) in the whole region 

G [1; 10] GeV^ with slightly higher values in the case of the improved Gaussian model. 

The SR result for the pion form factor is defined numerically as the average value of the 
right-hand side of SR flTUl) with respect to the Borel parameter g [Ml, M^]: 

1 r^^+ 

F^'^iQ') = M2_M2 J^^^ F{Q\M\som)dM\ (17) 

The obtained predictions for the pion form factor for both Gaussian NLC models with 
Ag = 0.4 GeV^, are shown in the right panel of Fig. [1] as dashed and solid curves, respec- 



tively, in comparison with the lattice result of [22j (dark grey strip limited from above at 
approximately 5 GeV^). These theoretical results are cornpared with the available experi- 



mental data 23, 24| and previous theoretical estimates [1, [g, 0]. We also show in this figure 
in the form of light grey strips the minimal theoretical uncertainties of the QCD SR results. 
The central curves of our predictions can be represented by the corresponding interpolation 
formulas: 

Fll,M' = ^ GeV^) = IMe-^-'^^'-'^x, (18a) 
^.^fmp(<3' = X GeV') = e-°-^2^"°'a; (O.OlGx^^ - 0.065a; + 0.58) , (18b) 
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FIG. 1: Left panel: Continuum threshold so(Q^) [GeV^] for the minimal (dashed line) and for the improved 
(solid line) NLC model. Right panel: Theoretical predictions for scaled pion form factor Q^F-iriQ^) obtained 
by different methods and models. Dashed line — minimal NLC model; solid line — improved NLC model (in 
both cases = 0.4 GeV^ has been used and the corresponding uncertainties are indicated by light gray 
strips delimited by similar lines). The followingdesignations are used: thick line between 1 and 4 GeV^ — 
standard QCD SR with local condensates dash-dotted line — LD QCD SR tri ang les — Cornell 

experimental data ,23]; diamonds — JLab experimental data 2J]. The recent lattice result 2^ is shown as 
a monopole fit containing error bars illustrated by a dark grey strip ending at « 4.5 GeV^. 



valid for ^ [i^ iq] GeV^, i. e., for x G [1, 10]. 



III. DISCUSSION AND CONCLUSIONS 



We calculated the electromagnetic pion FF using QCD SR with NLC jl, 1^ with a QCD 



vacuum nonlocality parameter = 0.4 GeV^ and using a perturbative spectral density 



1 



proposed in [l9|. This value receives support from a recent comprehensive analysis 
2^ of the CLEO data on the pion-photon transition.^ The use of NLC enables one to 



considerably enlarge the region of applicability of the QCD SR towards momenta as high 
as 10 GeV^ — in contrast to the standard QCD SR approach [sl, 0], where the SR can be 
applied only in the < 3 GeV^ region. 

The main conclusions of our investigation can be summarized as follows: 

• We found that the 0(Q;s)-contribution influences the pion form factor at the level of 
20%. This estimate is a little bit smaller than ones, obtained in (sl. [l9|. 

• We found that the central-line prediction of the improved model NLC model is inside 
the error strip of the minimal model up to = Q GeV^. Therefore, we may conclude 
that both models are equally good in this region. In view of the absence of more precise 
experimental data on the pion form factor at present, we cannot give any preference 
to one or the other of the two considered NLC models. 



^ Using somewhat higher values of this parameter, would entail a decrease of the pion form factor owing 
to a stronger influence of the nonlocality effects. The opposite effect appears for smaller values of this 
parameter and leads to an increase of the pion form factor. 
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It appears that our predictions are systematically higher than those of the LD ap- 
proach [5]. This can be easily understood in terms of the effective LD threshold 
So°(Q^). As we have already said, its value in the LD approach is well established 
only in the small-Q^ region. For higher values, it is not firmly fixed; for instance, in 
it was suggested to use a logarithmically increasing threshold 



l + a,(g2)/vr' 

which is ~ 0.62 GeV^ for ~ 3 GeV^. We estimated that in order to imitate the 
NLC QCD SR results in the LD approach, one needs to use Sq^{Q'^ = 10 GeV^) = 
0.87 GeV^. This means that the Sq^ uncertainty in the region of = 10 GeV^ is of 
the order of 30% This is the origin of the discussed difference between the LD results 
and ours. 

The lattice QCD results of 2^ are in excellent agreement with our predictions. 

Both, the minimal and the improved Gaussian model for the NLC give results which 
are in good agreement within errors with the currently available experimental data. 
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APPENDIX A: QCD SR PARAMETERS 



The parameters of the NLC are A = Xl/2 = 0.2 GeV^, {a,GG)/7v = 0.012 GeV^ and 
c^s {qq)^ = 1-83 ■ 10^^ GeV^. The nonlocal gluon-condensate contribution ^(^(M^) produces 
a very complicated expression. In analogy to the quark case, we model it by an exponential 



factor |8|,|JJ: $g(M' 



(M 



with 



0.4 GeV2. 

In order to determine the best value of the threshold Sq, we define the function for 
each value of and Sq as follows 



AT, 



M 



Nm 
1=0 



\i=0 



2n 



A^M + 1 



(A.l) 
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where we used = Ml +iAM, Am = {Ml - MI)/Nm, Nm = 10, and with e denoting 
the desired accuracy for — 1 (the actual value used in the computation is e = 0.07.) 
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